MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  Of  STANDARDS  1963-A 


Department  of  Ocean  Engineering 


CD 


Cambridge,  MA  02139 


’ ^ 
Q_ 
O 


Report  No.  78-5 

POTENTIAL  THEORY  OF  STEADY  MOTION  OF  SHIPS, 
PART  3:  WAVE  RESISTANCE 

by 

Francis  Noblesse 

November  1978 


D D C 


B 


This  work  was  supported  by  the  Office  of  Naval  Research, 
Contract  N00014-78-C-0169,  NR  062-525,  MIT  OSP  85949 


Approved  for  public  release;  distribution  unlimited. 


78  12  Od  . net  a 


SECURITY  CLASSIFICATION  OF  THIS  RACE  (Whmt  Dim  Kntmrmd) 


REPORT  DOCUMENTATION  PAGE 


12.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  NUMBER 


4,  TITLE  (and  Sabtttla) 


I.  TYPE  OP  REPORT  • PERIOO  COVERED 


POTENTIAL  THEORY  OF  STEADY  MOTION  OF  SHIPS , 
♦r  'PART  3#  WAVE  RESISTANCE^ 


^Technical  Report / 


«.  PERFORMING  ORO.  REPORT  NUMiER 


7.£UTHORflU 

"W  I 

Francis  INoblesse 


t.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Department  of  Ocean  Engineering  • 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


•I.  CONTROLLING  OFFICE  NAME  AND  AOORESS 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  S WORK  UNIT  NUMBERS 


Office  of  Naval  Research 


■ 

I 


November  1978 


TTOWBEROF  PAGES 

26 


MONITORING  AGENCY  NAME  » AOORESSflf  dlllarant  from  Controlling  Oltlea)  IS.  SECURITY  CLASS,  (at  this  ration) 

Unclassified 

TSiOECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


«.  DISTRIBUTION  STATEMENT  (ol  this  Raport) 


Distribution  of  this  document  is  unlimited 

— 


( ^7  P' 


DISTRIBUTION  STATEMENT  A, 

Approved  for  public  release] 
Distribution  Unlimited 


19.  KEY  WORDS  (Canthaaa  on  rararaa  slda  II  naaaaaarr  "h*  Idanllfr  »r  «•**  nrmtar) 

Ship  wave  resistance 
Wave  resistance  formula 
Ship  waves 
Potential  theory 


20.  ABSTRACT  (Canllnut  on  rararaa  alda  II  nacaaaaty  and  Idantltr  *F  numkaa) 

X 

s New  expressions  for  the  wave  resistance  of  a ship  in  steady  rectilinear 
motion  in  a calm  sea  are  presented.  V 


DO  ,ETn  1473  COITION  OF  I NOV  SI  IS  OBSOLETE 


S/N  010  2-014*  660 1 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  I 


POTENTIAL  THEORY  OF  STEADY  MOTION  OF  SHIPS, 
PART  3:  WAVE  RESISTANCE 


TABLE  OF  CONTENTS 


Page 


Abstract  i 

Introduction I 

1.  The  integral  equation  for  the  velocity  potential  <f>, 

and  the  approximations  and  <p^ 3 

2.  The  Kochin  free-wave  spectrum  function  fi(9) 

and  the  Havelock  wave-resistance  formula  7 

3.  Elementary  wave  resistance  formulas:  the  zeroth  approximation  Rq 

and  particular  limiting  cases  10 

4.  The  first  approximation  R^ 16 

5.  The  second  approximation  R£ 19 

6.  Discussion 21 

7.  Comparison  between  the  zeroth  approximation  R^  and  the  Michell 

approximation  R^  for  a wedge-like  ship-bow  form  24 

References 26 

Acknowledgments  26 


1 ACCESSION  (m 

ny 

NTIS  V.  ' 3 section  vn 

0DC  Cvii  Section  □ 

UNANNOUNCED 
JUSTIFICATION 

□ 

BY 

1 BISTRIBUTIOH/MAHABlllTK  COOK  J 

SFECIAL 

wz 

| 

i 


ABSTRACT 


New  expressions  for  the  wave  resistance  of  a ship  in  steady  rectilinear 
motion  in  a calm  sea  are  presented.  Specifically,  a sequence  of  iterative 
approximations  Rq,  R^,  R£,  ...  is  defined.  Of  particular  interest  are  the 
zeroth  approximation  Rq  and  the  first  approximation  R^ . The  zeroth  approximation 
Rq  provides  a new  simple  explicit  wave-resistance  formula  which  may  be  regarded 
as  a generalization  of  the  classical  formulas  proposed  by  Michell  in  1896  and 
Hogner  in  1932.  A noteworthy  feature  of  the  approximation  Rq  is  that  it  involves 
a line  integral  along  the  ship  waterline,  which  causes  a drastic  reduction  in  the 
wave  resistance  at  low  Froude  number  and  is  particularly  significant  for  blunt 
ship  forms.  Although  the  first  approximation  R^  is  of  course  more  complex  than 
the  ziioth  approximation  Rq,  it  provides  a fairly  simple  explicit  approximation 
to  the  wave  resistance  which  is  more  refined  than  the  approximation  Rq  and  may 
be  of  greater  usefulness  for  practical  purposes.  The  first  approximation  R^  is 
indeed  regarded  as  the  main  result  of  this  study.  Comparison  between  the  zeroth 
approximation  Rq  and  the  classical  Michell  approximation  for  a wedge-like  ship 
bow  form  suggests  that  the  present  theory  may  remedy  to  some  of  the  typical 
discrepancies  between  experimental  and  theoretical  wave-resistance  curves. 
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INTRODUCTION 


The  single  most  important  goal  of  a theory  of  steady  rectilinear  motion  of  a 
ship  in  a calm  sea  no  doubt  is  the  prediction  of  the  drag,  so-called  wave  resistance, 
experienced  by  a ship  as  a result  of  the  surface  gravity  waves  it  creates.  Exact 
expressions  (to  be  sure,  within  the  limitations  of  potential -flow  theory)  and  re- 
lated explicit  approximations  for  the  wave  resistance  of  a ship  are  presented  in 
this  study,  which  corresponds  to  Part  3 of  the  theory  of  steady  motion  of  ships 
developed  in  Noblesse  [1],  Only  the  case  of  displacement  ships  is  considered  here 
for  shortness;  however,  modification  of  the  analysis  and  of  the  final  formulas 
for  the  case  of  other  types  of  ships,  e.g.  fully-submerged  bodies,  multihull  vessels, 
and  surface-effect  ships,  is  straightforward  (problems  associated  with  lift  and 
cavitation  for  hydrofoils,  and  planing  effects  and  spray  formation  for  fast  boats, 
are  not  considered  in  the  present  theory  however) . 

The  wave  resistance  is  defined  by  a set  of  three  equations,  namely,  (i)  the 
classical  Havelock  wave-resistance  formula  (9)  expressing  the  wave  resistance  R in 
terms  of  the  Kochin  free-wave  spectrum  function  0(9),  which  is  directly  related  to  the 
free-wave  pattern  trailing  far  behind  the  ship,  (ii)  formula  (7)  defining  the  Kochin 
spectrum  function  0(9)  in  terms  of  the  velocity  potential  <j>  of  the  disturbance  flow 
caused  by  the  ship  in  its  "near  field",  and  (iii)  the  integral  equation  (1)  for  deter- 
mining the  velocity  potential  (j>  (this  integral  equation  was  derived  in  Part  2 of  [1]). 
Actually,  to  these  three  equations  one  should  add  equations  for  determining  the  hydro- 
dynamic  lift  and  moment,  and  the  resulting  sinkage  and  trim,  experienced  by  the  ship; 
however,  these  additional  equations  are  not  considered  explicitly  in  this  study. 

In  practice,  the  integral  equation  (1),  and  formulas  (7)  and  (9)  may  be  used  for 
defining  a sequence  of  iterative  approximations  to  the  wave  resistance.  Specifically, 


corresponding  to  the  successive  iterative  approximations  (|>q  = 0,  4>2*  •••  to  the 

solution  <J>  of  the  integral  equation  (1) , we  may  readily  associate  the  approximations 


(2q,  (2^,  &2»  •••  to  Kochin  free-wave  spectrum  function  (2(9),  and  the  approximations 
Rq,  R^,  Rj,  ...  to  the  wave  resistance  R,  by  using  formulas  (7)  and  (9)  in  which  we 


need  only  replace  <p  by  <p^>  (2  by  (2^,  and  R by  R^,  with  k = 0,  1,  2, 


the  iterative 


approximations  <J>j  and  <p^  to  the  solution  of  the  integral  equation  (1)  were  obtained 


previously  in  Part  2 of  [1],  and  are  given  in  this  Part  3 by  formulas  (2)  and  (3). 


The  approximations  (2^,  £2^,  and  (2 ^ are  derived  explicitly,  and  discussed  in  some 
detail.  The  general  recurrence  relations  defining  the  iterative  approximations  ^ and 
(2^  are  also  given  [specifically  by  formulas  (30)  and  (31)].  Of  main  interest  for 
practical  purposes,  however,  are  the  "zeroth  approximation"  (2^  and  the  first  approxi- 
mation (2^  corresponding  to  the  approximations  4>q  = 0 and  <J>j,  respectively. 
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The  zeroth  approximation  £2^(9),  given  by  formula  (10),  provides  a new  simple 
explicit  wave-resistance  formula,  which  may  be  regarded  as  a generalization  of  the 
classical  formulas  proposed  by  Michell  in  1898  and  Hogner  in  1932.  Indeed,  these 
classical  approximations  are  obtained  in  the  present  study  as  particular  limiting 
cases  of  the  zeroth  approximation  £2q.  A noteworthy  feature  of  the  approximation  £2^ 
is  that  it  involves  a line  integral  along  the  ship  waterline,  as  it  may  readily  be 
seen  from  formula  (10).  This  line  integral  is  shown  to  cause  a drastic  reduction  in 
the  value  of  the  wave  resistance  at  low  Froude  number,  and  is  particularly  significant 
for  blunt  ship  forms. 

The  first  approximation  £2^(9),  defined  by  formula  (22),  is  associated  with  the 
approximation  given  by  formula  (2).  Although  this  first  approximation  £2^  is  ob- 
viously more  complex  than  the  zeroth  approximation  £2^,  it  provides  a fairly  simple 
explicit  approximation  to  the  Kochin  spectrum  function  (2(9)  which  clearly  is  more 
refined  than  the  approximation  £2q,  and  may  be  of  greater  usefulness  for  practical 
purposes.  In  particular,  the  first  approximation  £2^  incorporates  correction  terms 
associated  with  the  hull  boundary  condition,  including  effects  of  sinkage  and  trim, 
and  corrections  for  free-surface  nonlinearities.  For  practical  purposes,  the  first 
approximation  is  regarded  as  the  main  result  of  the  present  study.  Various  ways  of 
improving,  and  simplifying,  this  approximation  are  briefly  discussed. 

Finally,  preliminary  numerical  results  are  presented.  Specifically,  a comparison 
between  the  zeroth  approximation  Rq  and  the  classical  Michell  approximation  for  a 
wedge-like  ship  bow  form  is  made.  The  results,  shown  in  figure  2,  indicate  that  — 
compared  with  the  Michell  approximation  (M)  — the  new  approximation  yields  a 
significant  overall  reduction  in  the  value  of  the  wave  resistance,  a very  appreciable 
reduction  in  the  magnitude  of  the  oscillations  (humps  and  hollows)  in  the  wave 
resistance  curve,  and  a notable  phase  shift  of  these  oscillations  towards  lower 
values  of  the  Froude  number.  These  encouraging  preliminary  results  suggest  that 
the  present  theory  may  thus  remedy  to  some  of  the  typical  discrepancies  between  ex- 
perimental and  theoretical  wave  resistance  curves. 
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1.  The  integral  equation  for  the  velocity  potential  <j),  and  the  approximations  4^  and  ^ 


We  consider  a displacement  ship  in  steady,  rectilinear  motion  at  the  free 
surface  of  an  otherwise  calm  sea  assumed  to  be  of  infinite  depth  and  lateral 
extent.  Water  is  supposed  to  be  homogeneous  and  incompressible.  Surface  tension 
is  neglected.  Irrotational  flow  is  assumed;  however,  ad  hoc  corrections  for  the 
viscous  boundary  layer  and  wake  around  and  behind  the  ship  are  incorporated  into 
the  present  potential-flow  model,  in  the  manner  explained  in  [1].  Ad  hoc  cor- 
rections for  effects  of  spray  formation  along  the  waterline,  and  for  effects  of 
wavebreaking  are  also  included,  as  in  [1]. 

A moving  system  of  coordinates  attached  to  the  ship  is  chosen,  so  that  the  flow 
is  independent  of  time.  The  z axis  is  taken  vertical,  positive  upwards,  with  the 
undisturbed  free  surface  taken  as  the  plane  z = 0.  The  x axis  is  parallel  to  the 
direction  of  motion  of  the  ship  and  positive  towards  the  ship  stern.  Flow  variables 
are  rendered  dimensionless  with  respect  to  the  speed  of  the  ship  U,  the  acceleration 
of  gravity  g,  and  the  fluid  density  p,  as  it  is  shown  in  equations  (1.5)  in  [1]. 

The  hydrodynamical  problem  amounts  to  determining  the  (dimensionless)  velocity 
potential  <p  of  the  disturbance  flow  caused  by  the  ship.  This  problem  was  formulated 
in  [1]  as  a "generalized  Neumann-Kelvin  problem"  in  a "solution  domain"  (d)  bounded 
by  some  "fictitious  hull  surface"  (h) , which  may  (but  need  not)  be  taken  as  the 
submerged  hull  of  the  ship  in  position  of  rest  [in  principle,  the  "fictitious  hull 
surface"  (h)  may  be  chosen  arbitrarily] , and  the  undisturbed  free  surface  (f ) , which 
is  the  portion  of  the  plane  z = 0 outside  the  intersection  curve  (c)  of  the  surface 
(h)  with  the  plane  z = 0.  It  is  shown  in  [1]  that  the  "generalized  Neumann-Kelvin 
problem"  can  be  formulated  in  "integral  form",  given  by  the  integral  equation  (2.21). 
This  integral  equation  is  given  here  for  easy  reference: 


<J>0  = j Gvda  -“*1)  Gv2yds  + J (4>-<f>0)Gnda  + *V  [G(a4>s+r<J>t)-(<}>-<J>0)Gx]yds 
> h ' c ' h c 


+ j G(qnl+qwb)dxdy  + J^q^+q^da  - ^G[  (q^^vy-q^ds  + J ^dv 


' 

- 1 


't 


(1) 


where  the  significance  of  the  various  undefined  symbols  will  now  be  explained. 

The  symbol  4>q  is  meant  for  4>(Xq)  where  Xq  is  an  arbitrary  point  in  the  solu- 
tion domain  (d)  including  its  boundary  (h)  + (f)  + (c) , while  4>  is  meant  for  4>(x) 
where  x represents  the  "point  of  integration"  (integration  variable)  in  the  above 
integrals;  the  point  thus  is  the  "field  point"  where  the  potential  is  being 
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evaluated,  while  x represents  the  "dummy"  variable  of  integration.  The  function 
G = G(Xq,x)  is  the  fundamental  solution  (Green  function)  appropriate  for  the  problem. 
Specifically,  the  function  G(Xq,x)  represents  the  (dimensionless)  linearized  velocity 
potential  of  the  disturbance  flow  caused  at  point  x^(x^,y^,ZQ<0)  by  a unit  "outflow" 
at  point  x(x,y,z<0),  associated  with  a submerged  source  if  z < 0 or  a flux  across 
the  free  surface  if  z = 0,  in  an  oncoming  uniform  stream  along  the  positive  x axis. 
Simplified  new  expressions  for  the  fundamental  function  G and  its  gradient  VG  may 
be  found  in  Noblesse  [2], 

The  symbol  da  in  the  three  surface  integrals  over  (h)  represents  the  differ- 
ential element  of  area  of  (h) , while  ds  in  the  three  line  integrals  around  the 
"waterline"  (c)  represents  the  differential  element  of  arc  length  of  (c) , and  dv 
in  the  volume  integral  over  the  solution  domain  (d)  is  the  differential  element  of 
volume.  In  the  first  integral,  we  have  v = v(x)  = n(x)*i,  where  n(x)  is  the  unit 
inward  (that  is,  n is  pointing  towards  the  interior  of  the  ship)  normal  vector  to  (h) 
at  point  x of  (h) , and  i is  the  unit  positive  vector  along  the  x axis.  In  the  line 
integrals  around  (c),  we  have  v = v(s)  5 n(s) *i  where  n(s)  is  the  normal  to  (h)  at 
point  s of  (c) , while  p is  defined  as  p = p(s)  = n'(s)‘t,  where  n'  is  the  unit  inward 
normal  vector  to  (c)  in  the  plane  z =»  0;  in  the  (fairly  common)  case  when  the  sur- 
face (h)  intersects  the  plane  z =*  0 orthogonally,  we  have  n(s)  = n'(s)  and  v(s)  = 

p(s).  In  the  second  line  integral  around  (c),  the  symbols  a and  T are  defined  as 
->  “►  ■+■  ->  . 

a = s*i  and  x = t*i,  where  s = s(s)  is  the  unit  tangent  vector,  at  point  s,  to  the 
"waterline"  (c)  oriented  in  the  counterclockwise  direction  in  the  (x,y)  plane,  and 
t = t(s)  is  the  unit  vector  tangent  to  (h) , mutually  orthogonal  to  s(s)  and  the 
normal  n(s)  to  (h)  at  point  s of  (c) , and  pointing  downwards.  In  this  line  integral, 

3<j>(s,t,n)/3t,  and  G 


the  notation  <f>s  = 3<t> (s, t,n) /3s, 


3G(Xq,x)/3x  was  used  for 


shortness.  The  usual  notation  G^  = VG(Xq,x> *n(x)  was  also  used  in  the  third  integral. 
It  will  be  noted  that  the  axes  x,  y,  and  z,  the  "fictitious  hull  surface"  (h) , the 
undisturbed  free  surface  (f ) , the  "waterline"  (c) , the  elements  of  area  da  and  of 
arc  length  ds,  and  the  unit  vectors,  i,  n,  n',  and  i:,  are  shown  in  figure  1. 

It  remains  to  define  the  terms  q^,  q^f,  qg,  qfel,  q^,  and  q^.  The  term  qQ^ 
in  the  integral  over  the  undisturbed  free  surface  (f)  represents  the  "NonLinear 
free-surface  correction  flux",  which  accounts  for  the  nonlinear  terms  in  the  free- 
surface  boundary  condition  and  for  the  difference  in  position  between  the  actual 
free  surface  and  the  plane  z = 0 of  the  undisturbed  free  surface  where  the  free- 
surface  condition  is  enforced  for  mathematical  simplicity;  the  "nonlinear  free- 
surface  flux"  q^  is  given  by 
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Hnl 


<P  + <t> 


XX 


ClV4>|2)' 


- <j>  - i-i 
yx  2 1 


xx 


(la) 


z = 0 


The  term  q is  meant  for  the  "Hull-Form  correction  flux"  associated  with  the  fact 
that  the  "fictitious  hull  surface"  (h)  may  be  different  from  the  actual  ship  hull 
surface  (H) ; the  "hull-form  correction  flux"  is  given  by 


qhf  = (i  + V<p)H  • N - (i  + V<J>)h  • n , 


(lb) 


where  N is  the  unit  inward  normal  vector  to  (H) , and  the  notation  ( )u  and  ( ),  means 

n h 

that  the  expression  between  the  parentheses,  namely  V<j>,  is  to  be  evaluated  on  (H) 
and  (h) , respectively.  Finally,  the  terms  q^p  qg,  and  q^  are  ad  hoc  corrections 

associated  with  the  viscous  Boundary  Layer  around  the  ship  hull.  Spray  formation 
along  the  waterline,  WaveBreaking  at  the  free  surface,  and  the  viscous  Wake  trailing 
behind  the  ship,  respectively. 

It  may  be  useful  to  emphasize  that  the  integral  equation  (1)  is  valid  for  Xq 
in  the  solution  domain  (d)  including  its  boundary  (h)  + (f)  + (c),  as  it  was  already 
noted,  so  that  this  equation  defines  the  disturbance  velocity  potential  <p  everywhere 
in  (d)  + (h)  + (f)  + (c) . The  solution  <J>  of  the  integral  equation  (1)  may  be  deter- 
mined in  practice  by  using  an  iterative  method  of  solution,  as  it  is  explained  in  [1]. 

Briefly,  if  we  merely  ignore  all  the  unknown  terms  in  the  integral  equation  (1), 
we  obtain  the  initial  approximation  <J>^(Xq),  which  is  thus  given  by 


(2) 


A second  approximation  can  then  be  determined  by  evaluating  the  previously- 

ignored  unknown  terms  (however,  the  correction  terms  for  "real-fluid  effects"  q^p 
qg,  q^*  and  qw  are  neglected  in  the  following  second  approximation  <t>^)  on  the  basis 
of  the  Initial  approximation  <j>j.  to  i}>.  After  some  manipulations,  we  may  obtain 


(3) 


where  the  notation  (v-hf)1)  is  meant  for  (v)  +(4>^)„  = i*N+(V4>  ) *N  = (±+V<J>_)  -N  > 

nil  Hnn  in 

which  represents  the  fluid  flux  across  the  actual  ship  hull  (H)  in  the  initial 
approximation  i + V<J>j,  and  the  nonlinear  free-surface  flux  qR^  is  given  by 


p 


6 


■ 


v-  n^i2)*  -♦X  + O, 


(3a) 


While  higher  iterative  approximations  <f>^ , k ^ 3,  could  readily  be  defined  in 
principle,  numerical  evaluation  of  these  higher  approximations  would  probably  be 
too  considerable  a task  in  practice.  However,  it  seems  reasonable  to  hope  that 
the  second  approximation  (J^,  or  even  perhaps  the  initial  approximation  may  be 
sufficiently  accurate  for  most  practical  applications.  Further  discussion  of  the 
above  approximations  and  may  be  found  in  [1],  where  variations  about  these 
iterative  approximations  are  also  examined.  : may  finally  be  noted  here  that  the 
value  of  the  potential  (J)  in  any  portion  of  the  fluid  domain  (D)  outside  the  solution 
domain  (d) , e.g.  in  the  region  between  the  actual  free  surface  (F)  and  the  plane 
z = 0 wherever  (F)  lies  above  the  plane  z = 0,  may  be  determined  from  the  value  of 
<t>  in  (d)  by  means  of  analytical  continuation. 


Figure  1:  Definition  sketch 
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2.  The  Kochin  free-wave  spectrum  function  ft(9)  and  the  Havelock  wave-resistance  formula 


A ship's  wave  resistance  is  related  in  a simple  manner  to  the  wave  pattern 
trailing  behind  the  ship,  as  it  is  well  known  since  Havelock  [3].  The  disturbance 
velocity  potential  <j>  in  the  "far  field",  that  is  at  a large  distance  away  from  the 
ship,  may  be  expressed  in  terms  of  the  value  of  <|>  in  the  "near  field",  that  is  on 
the  hull  surface  (h)  + (c)  and  on  the  free  surface  (f)  in  the  vicinity  of  the  ship, 
by  means  of  the  equation 


<J>Q  - Gvda  - Gv^yds  + | <}>G^da  + [G  (0<}>g+T<J>t) -<{>Gx]yds  + G(q  b+q^b)dxdy  + 

+ G(qhf+qbl)da  - ^ G[  (qhf+qbl)vy-qg]ds  + Gq^dv  » (4) 


which  may  be  derived  from  the  integral  equation  (1)  by  replacing  the  term  <p  ~ <Pq 
in  the  third  and  fourth  integrals  by  <j)  since  we  have  ^ 0 as  |xq|  -*■  °°  , and 

therefore  |<(>q|  « |<f>|  for  Xq  in  the  far  field  and  x in  the  near  field.  It  is 
interesting  to  note  here  in  passing  that  while  expression  (4)  for  the  disturbance 
velocity  potential  4>(Xq''  in  the  far  field  was  obtained  above  as  the  "far-field 
limit"  of  the  integral  equation  (1) , it  is  actually  valid  also  for  Xq  in  the  near 
field,  provided  only  that  Xq  is  strictly  outside  the  hull  surface  (h)  + (c) , that 
is  equation  (4)  holds  for  Xq  in  the  domain  (d)  + (f)  - (h)  - (c) ; indeed,  equation 
(4)  corresponds  to  equation  (2.11a)  in  [1]. 

Far  behind  the  ship,  that  is  for  Xq  » 1,  the  foregoing  expression  for  <p^ 
can  actually  be  greatly  simplified  by  replacing  the  Green  function  G(Xq,x)  by  the 
well-known  asymptotic  approximation 


G ^ Re 


{ (zQ+z)+i[ (xQ-x)cos0+(yo-y)sin0 ] }sec  0 2 


sec  0d0  as  x^  -*■  00 


which  may  be  obtained,  for  instance,  from  equations  (19),  (23)  and  (10c)  in 
Noblesse  [4];  the  above  approximation  differs  from  the  exact  expression  for  G by 
terms  associated  with  a near-field  disturbance  which  "dies  out"  like  1/Xq  as 
Xq  -*■  00  , and  may  be  discarded  inasmuch  as  we  are  only  interested  in  the  trailing 
wave  pattern  far  behind  the  ship. 


8 


By  using  the  asymptotic  approximation  (5)  into  equation  (4) , we  may  express 
the  disturbance  potential  <J>(Xq)  far  behind  the  ship  in  the  form 


4>(x  ) ^ Re  ^ JI(0)sec  9 e 


2 [z^+i (x^cosQ+y^sinS) ]sec  9 


d0  as  xn  -*■  °°  , (6) 


where  the  function  £2(9)  is  given  by 


[E°(CTc(>s+T<{>t)-<j>E°]yds  + 


£2(9)  = Evda  - E^v^yds  + (fiE^da  + *6  [E^(ct<J>  +tc{)  )-(j)E^]yds  + 

'h  ' c -* h ) c s x 

E°(qnl+qwb)dxdy  + | E(qhf+qbl)da  “f  E°t^hf+qbl)VlJ“qs]ds  + j/V*' 
' h J c J a 


, (7) 


with  the  functions  E = E(9;x)  and  EU  = E(0;x,y,z=O)  defined  as 


E = e 


[z-i(xcos9+ysin0) ]sec  9 


fc  = e 


-i(xcos9+ysin0)sec  9 


(7a, b) 


The  equation  of  the  free  surface  far  behind  the  ship  is  given  by  Zq  = -<J)x(xQ,yg,ZQ=0) 

since  the  nonlinear  terms  in  the  free-surface  boundary  condition  may  be  neglected 
at  a sufficiently  large  distance  away  from  the  ship.  By  using  equation  (6)  we  may 
then  obtain 


zn  * Re  f 
0 TT 


£2(9)sec  0 e 


^ i(xQCOs0+ygSin0)sec  0 


d9  as  Xq  -*■  00  . 


Equations  (6)  and  (8)  express  the  potential  (J)(xq)  and  the  equation  of  the  free 
surface  ZQ^.y^)  far  downstream  from  the  ship  in  terms  of  a familiar  superposition 
of  elementary  plane  waves  with  amplitude  £2(9)sec^0.  The  function  £2(9)  corresponds 
to  a particular  case  of  the  well-known  function  introduced  by  Kochin  for  determining 
the  drag,  lift,  and  moment  acting  upon  a ship,  which  is  often  referred  to  as  the 
"free-wave  spectrum"  in  the  literature  on  "wave  analysis"  (see  for  instance  Eggers, 
Sharma,  and  Ward  [5]);  the  function  £2(9)  defined  by  equation  (7,a,b)  will  be  referred 

to  as  the  "Kochin  free-wave  spectrum  function"  in  the  present  study. 

2 6 

The  (dimensionless)  wave  resistance,  R say  (R  = R'g  /pU  , where  R'  is  dimen- 
sional), may  be  directly  determined  from  the  "Kochin  free-wave  spectrum  function" 

£2(9)  by  means  of  the  well-known  "Havelock  wave -resistance  formula" 


9 


if  'K 

R = it-  |(2(9)  | secJ9de 


which  is  given  for  instance  in  [5]  equation  (12)  p.  118.  In  the  usual  case  of  a 
ship  which  is  symmetric  with  respect  to  its  centerplane  y = 0,  the  function  |(2(0)| 
is  even,  so  that  the  "Havelock  wave-resistance  formula"  (9)  may  be  expressed  in  the 


R-i 

TT 


|n(9) |2sec39d9 


The  "Havelock  wave  resistance  formula"  (9)  for  the  wave  resistance  R,  formula 
(7)  for  the  "Kochin  free-wave  spectrum  function"  0(9),  and  the  integral  equation  (1) 
for  determining  the  disturbance  velocity  potential  4>  in  the  "near  field"  of  the  ship, 
form  a set  of  three  equations  for  determining  the  wave  resistance  of  a ship  in 
steady  motion.  Actually,  to  these  three  equations  we  must  also  add  equations  for 
determining  the  sinkage  and  trim  experienced  by  the  ship,  which  requires  evaluation 
of  the  hydrodynamic  lift  and  moment  acting  upon  the  ship.  However,  the  equations 
for  lift  and  moment  and  sinkage  and  trim  will  not  be  considered  explicitly  in  the 
present  study,  which  is  mainly  concerned  with  the  "hydrodynamical  problem"  defined 
in  Part  1 of  [1],  that  is  the  problem  of  predicting  the  flow  caused  by  a ship,  and 
in  this  Part  3 the  wave  resistance,  assuming  the  shape,  and  position,  of  the  ship 
to  be  given. 

The  integral  equation  (1)  and  formulas  (7)  and  (9)  are  essentially  "exact", 
that  is  within  the  limitations  of  the  present  potential-flow  model,  and  in  principle 
these  equations  thus  provide  a basis  for  obtaining  an  essentially  "exact"  value  of 
the  wave  resistance  of  a ship  in  steady  motion.  In  practice,  equations  (1),  (7), 
and  (9)  may  be  used  for  determining  a sequence  of  "iterative  approximations"  to  the 
wave  resistance  R:  specifically,  corresponding  to  the  successive  iterative  approxi- 
mations 4>q  = 0,  <j)^,  4> ^ » •••  to  the  solution  <p  of  the  integral  equation  (1),  we  may 
readily  associate  the  approximations  (2^,  (2^,  •••  to  the  Kochin  free-wave 

spectrum  function  (2(0),  and  the  approximations  Rq,  R^,  R^,  ...  to  the  wave  resistance 
R,  by  using  formulas  (7)  and  (9)  in  which  we  need  only  replace  <j>  by  (j)^,  (2  by  (2^,  and 
R by  R^,  with  k = 0,  1,  2,  ...  . The  approximations  (2^,  (2^,  and  ^2  are  given  and 
discussed  in  some  detail  in  the  following  three  sections;  several  particular  limiting 
cases  of  these  approximations  are  also  discussed. 
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3.  Elementary  wave-resistance  formulas:  the  zeroth  approximation  R^  and  particular 
limiting  eases 


In  this  section,  five  elementary  wave  resistance  approximations  are  derived 
and  discussed.  These  are  the  "zeroth  approximation"  Rq,  and  its  four  particular 
limiting  cases  the  Hogner  "fine-ship  approximation"  R^,  the  Hogner  "flat-ship 
approximation"  R^,  the  Michell  "thin-ship  approximation"  R^,  and  the  Maruo-Tuck- 
Vossers  "slender-ship  approximation"  R^^. 

The  simplest,  and  no  doubt  the  crudest,  possible  approximation  to  the  "Kochin 
free-wave  spectrum  function"  £2(9),  and  hence  to  the  wave  resistance  R,  is  obtained 
if  one  merely  ignores  the  integral  equation  (1)  for  determining  the  disturbance 
velocity  potential  <J>  in  the  "near  field",  and  puts  <p  = 0 in  formula  (7)  for  £2(9). 
This  approximation  may  perhaps  be  regarded  as  a "slender-ship  approximation"  since 
it  can  be  expected  to  be.  the  better  the  smaller  the  disturbance  potential  <p , that 
is  the  more  "slender"  the  ship;  however,  it  will  simply  be  referred  to  as  the 
"zeroth  approximation"  in  the  present  study.  This  "zeroth  approximation"  thus  is 
given  by 


v»  ■ L 


[z  - i(xcos9  + ysin0)]sec  9 


v(x)da(x) 


n 

-i 


-i(xcos9  + ysin9)sec  9 ^ 


V (s)u(s)ds 


(10) 


as  it  may  readily  be  obtained  from  formulas  (7,a,b) . 

If  the  ship  is  sufficiently  "fine",  that  is  if  the  angle  between  the  "water- 

2 

line"  (c)  and  the  x axis  is  sufficiently  small,  we  have  V |y|«|v|«l  and  the  line 
integral  in  formula  (10)  may  be  neglected  in  comparison  with  the  surface  integral. 
It  is  interesting  that  the  approximation,  £2^(9)  say,  defined  by  this  surface 
integral,  that  is 


£2h(9) 


[z  - i(xcos9  + ysin9)]sec  9 


v(x)da(x) 


(ID 


in  fact  corresponds  to  the  approximation  proposed  by  Hogner  [6]  in  19j2;  indeed, 
the  wave  resistance  approximation  R^  obtained  by  replacing  £2(9)  by  £2^(9)  in  formula 
(9)  is  identical  to  Hogner's  wave  resistance  formula.  The  approximation  £2^(9),  R^ 
will  be  referred  to  as  the  Hogner  "fine-ship  approximation",  or  simply  as  the 
Hogner  approximation. 
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If  the  equation  of  the  (fictitious)  hull  surface  (h)  may  be  written  in  the 

2 

form  y = ± b(x,z),  where  b(=  gB/U  ) is  dimensionless  like  the  other  flow  variables, 
and  the  points  (x,z)  belong  to  the  projection  (h^)  of  the  surface  (h)  onto  the  ship 
centerplane  y = 0,  the  surface  integral  in  formula  (11)  may  be  transformed  into  the 
following  double  integral 


tty(0) 


l 


zsec  9 - ixsec8 


cos[b(x,z)tan9sec9]bx(x,z)dxdz 


(11a) 


If  the  ship  is  "thin",  that  is  if  b(x,z)  is  sufficiently  small  that  the  term 
cos[b(x,z)tan0sec9]  may  be  approximated  by  1,  the  Hogner  "fine-ship  approximation" 
tty(0)  given  in  equation  (11a)  becomes  the  famous  Michell  "thin-ship  approximation" 
5^(0)  say,  first  obtained  by  Michell  in  1898  and  given  by 


V9) 


zsec  0 - ixsec9 


bx(x,z)dxdz 


(12) 


Several  differences  between  the  Hogner  approximation  tty,  Ry  and  the  Michell 

approximation  RM  are  readily  apparent  and  may  be  noted  here.  An  interesting 

difference  between  tty(0)  and  tt^(9)  is  that  while  the  Michell  approximation  0^(9)  is 

proportional  to  the  ship's  beam  B,  and  R^  therefore  is  proportional  to  B^  (as  it  is 

very  well  known) , the  Hogner  approximation  tty(0)  clearly  is  not  proportional  to  B 

(although  tty  becomes  proportional  to  B in  the  limit  B -*■  0)  . It  can  also  readily  be 

seen  that  differences  between  Ry  and  R^  may  be  expected  to  be  the  larger  the  bigger 

the  beam  and  the  smaller  the  Froude  number;  this  is  due  to  the  fact  that  we  have 

b = gB/U^  = (B/L)(gL/U^)  = e/F^,  where  e = B/L  is  the  "geometrical  thinness"  ratio, 
1/2 

and  F = U/(gL)  is  the  Froude  number  based  on  the  ship  length  L.  It  is  interesting 

to  note  that  the  "thin-ship  assumption"  b«l  used  in  deriving  expression  (12)  for 

the  Michell  approximation  tt^  from  expression  (11a)  for  the  Hogner  approximation  tty 

implies  not  only  "geometrical  thinness",  characterized  by  e = B/L«l  , but  also 

2 

"Froude  thinness",  that  is  e/F  «1.  It  may  also  be  interesting  to  note  that  while 
most  ships  would  seem  to  be  "geometrically  thin",  in  that  e for  most  ships  is  fairly 
small  compared  with  1,  they  are  usually  not  "Froude  thin",  in  that  the  value  of  the 
ratio  e/F  is  commonly  equal  to  1 or  2,  which  hardly  seems  small  enough  to  justify 
approximating  the  term  cos[b(x,z)tan0sec9]  by  1.  However,  this  approximation  is 
certainly  justified  for  small  values  of  9,  corresponding  to  the  "transverse  waves" 
in  the  free-wave  spectrum;  on  the  other  hand,  the  approximation  btan9sec0«l  clearly 
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is  not  valid  for  values  of  | 0 | close  to  tt/2,  that  is  for  the  "divergent  waves"  in 
the  wave  spectrum.  Two  factors  which  contribute  to  reducing  the  differences  between 
and  may  finally  be  noted:  these  factors  are  that  the  term  b(x,z)tan0sec9  is 
small  for  b(x,z)  small,  that  is  at  the  ship  bow  and  stern  (which,  however,  are  the 
main  contributors  to  the  wave  resistance),  and  for  9 small,  that  is  for  the  part  of 
the  wave  resistance  corresponding  to  the  "transverse  waves"  (which  are  known,  how- 
ever, to  account  for  a significant  portion  of  the  wave  resistance). 

If  the  equation  of  the  (fictitious)  hull  surface  (h)  is  expressed  in  the  form 

2 

z * - d(x,y),  where  d(=gD/U  ) is  dimensionless,  and  the  points  (x,y)  belong  to  the 
projection,  (hz)  say,  of  the  surface  (h)  onto  the  "waterplane"  z = 0,  the  surface 
integral  (11)  may  be  transformed  into  the  following  double  integral 


flgO) 


~[d(x,y)  + i(xcos9  + ysin9)]sec  9 


d (x,y)dxdy 


If  the  ship  is  "flat",  that  is  if  d(x,y)  is  sufficiently  small  that  the  term 
2 

exp[-d(x,v)sec  9]  may  be  approximated  by  1,  we  obtain  the  Hogner  "flat -ship 
approximation",  £2^(9)  say,  given  by 


i4<9)  ■ IL 


-i(xcos9  + ysin9)sec  9 


dx(x,y)dxdy 


(13) 


The  approximation  to  the  wave  resistance  associated  with  the  above  approximation  . 
£2^  to  the  Kochin  free-wave  spectrum  function  is  identical  to  the  "flat-ship  wave 
resistance  formula"  proposed  by  Hogner  [6]  in  1932.  The  differences  between  the 
Hogner  "fine-ship  approximation"  £2^  and  the  Michell  "thin-ship  approximation" 
discussed  previously  may  also  be  observed  — with  evident  modifications  — between 
the  Hogner  approximations  £2^  and  £2^  , In  particular,  it  will  be  noted  that  £2^  is 
proportional  to  the  draft  D of  the  ship,  and  the  wave  resistance  R^  therefore  is 
proportional  to  D“",  which  is  not  true  for  the  approximation  £2^,  R^.  It  will  also 
be  noted  that  the  approximation  d(x,y)sec^9«l  must  obviously  break  down  as  1 9 1 - tt/2 
and  as  the  Froude  number  F -*■  0. 

If  the  ship  is  both  "thin"  and  "flat",  that  is  if  the  term  exp (-iytan9sec9) 

f 2 

in  formula  (13)  for  the  Hogner  "flat-ship  approximation"  £2^  and  the  term  exp(zsec  9) 
in  formula  (12)  for  the  Michell  "thin-ship  approximation"  £2^  can  be  approximated  by 
1,  both  approximations  (12)  and  (13)  yield  the  Maruo-Tuck-Vossers  "slender-ship 
approximation" 


r 


13 


stern 


fiMTV(0) 


-ixsec0 


A(x)dx 


(14) 


bow 

where  A(x)  is  given  by 
+b(x) 


A(x) 


dx(x,y)dy 


b^(x,z)dz 


(14a) 


-b(x)  -d(x) 

with  d(x)  and  b(x)  representing  the  local  (at  section  x)  draft  and  half  beam  of  the 
ship,  respectively.  The  above  approximation  S^y  was  obtained  by  Maruo  [7],  Tuck  [8], 


and  Vossers  [9]  in  the  early  1960's  by  using  the  method  of  "matched  asymptotic 

expansions".  It  may  readily  be  seen  that  the  Maruo-Tuck-Vossers  approximation 

is  proportional  to  the  product  BD  of  the  ship's  beam  and  draft,  and  the  wave  resistance 

2 2 

is  proportional  to  B D . An  appealing  feature  of  the  approximation  R^y  defined 
by  formulas  (14, a)  and  (9)  resides  in  its  remarkable  simplicity.  Unfortunately,  this 
wave  resistance  approximation  is  practically  useless  for  realistic  ship  hull  forms. 

In  particular,  the  source  strength  A(x)  must  vanish  at  the  ship  bow  and  stern  for 
the  wave  resistance  integral  (9)  merely  to  exist;  this  "existence  condition"  is 


extremely  restrictive,  and  indeed  is  not  satisfied  by  usual  ship  hull  forms. 

It  may  be  interesting  to  note  that  the  approximations  R^,  R^,  and  R^^,  which 
were  obtained  in  the  foregoing  as  particular  limiting  cases  of  the  Hogner  approximation 


Rjj,  can  also  be  obtained  as  the  first-order  (linearized)  approximations  in  asymptotic 
expansions  corresponding  to  the  assumptions  that  the  ship  is  "thin",  "flat",  and 
"slender",  respectively,  by  performing  systematic  perturbation  analyses  starting 
from  the  usual  differential  formulation  of  the  problem  of  steady  motion  of  a ship, 
as  it  is  well  known.  In  other  words,  the  "perturbation  approximations"  R^,  R^,  and 
R^py,  are  embodied  into,  indeed  are  particular  limiting  cases  of,  the  Hogner  approxi- 
mation Rjj,  which  itself  corresponds  to  the  "fine-ship  limit"  of  the  "zeroth  approxi- 
mation" Rq  obtained  in  the  present  study. 

While  the  neglect  of  the  line  integral  in  formula  (10)  for  the  zeroth  approxi- 
mation 0q(9)  can  be  justified  in  the  case  of  a "fine"  ship,  there  is  no  a priori 
reason  for  neglecting  this  line  integral  in  the  case  of  a ship  with  a blunt  bow  or/and 
stern.  The  line  integral  is  also  particularly  significant  at  low  values  of  the  Froude 
number;  specifically,  the  line  integral  causes  a drastic  reduction  in  the  value  of  the 
wave  resistance  at  low  Froude  number,  as  it  will  now  be  shown.  For  this  purpose,  we 
begin  by  considering  the  surface  Integral  in  formula  (10),  that  is  the  Hogner  "fine- 
ship  approximation"  &g(®)  given  by  formula  (11)  . It  is  convenient  here  to  introduce 
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the  coordinates  £(5,0,0  which  are  rendered  dimensionless  in  terms  of  some  character- 
istic dimension  of  the  ship,  say  its  length  L,  i.e.  we  have  t = X/L  where  X is  dimen- 
sional.  By  using  the  relation  x 
number,  into  equation  (11),  we  obtain 


9 -►2  1 12 

Xg/U  = £/F  , where  F = U/(gL)  is  the  Froude 


nH(9) 


-4 


F - i(£cos0  + nsin9)]sec^9 


v(C)da(C) 


(15) 


= r.4 


where  da  = F da  is  the  differential  element  of  area  of  the  surface  (h)  in  terms  of 
the  dimensionless  variables  X, 

The  major  contribution,  as  F + 0,  to  the  surface  integral  (15)  stems  from  the 


immediate  vicinity  of  the  "waterline"  (c) , due  to  the  rapid  exponential  decay  of  the 


-2  2 

factor  exp(F  £sec  9)  as  F -*■  0.  For  simplicity,  we  will  restrict  our  attention  to 


immediate  neighborhood  of  the  waterline  (c)  may  then  be  expressed  in  the  form  da 

_ „2 


% p 


-4 


-2  2 
F ^sec  0 


d£  e 


i 


-2  2 
-iF  (£cos9  + nsin9)sec  9 


v (a)  da 


where  the  undefined  lower  limit  of  integration  in  the  first  integral  may  actually 

-2  2 

be  taken  as  (-°°)  since  the  term  exp(F  £.sec  9)  is  "exponentially  small"  as  F -►  0. 


We  may  finally  obtain 


ftg(9)  'v-  F ^ cos^9  e 


1 


-iF  ^(£cos9  + nsin9)sec^9 


U(a)da  as  F 


(16) 


"zeroth  approximation"  (1q(9) 


the  particular,  but  fairly  common,  case  when  the  (fictitious)  hull  surface  (h)  inter- 
sects the  plane  z = 0 orthogonally,  that  is  when  (h)  is  vertical  sided  in  the  vicinity 
of  the  waterline.  The  differential  element  of  area  da  of  the  surface  (h)  in  the 


dad;;,  where  da  = F ds  is  the  differential  element  of  arc  length  along  (c) . In  the 
low-Froude-number  limit,  the  surface  integral  (15)  may  then  be  approximated  as  follows 

0 


since  we  have  v(a)  = p(a)  in  the  present  case  when  (h)  intersects  the  plane  z = 0 
orthogonally. 

By  combining  the  line  integral  in  formula  (10)  and  the  line  integral  (16) , which 
is  a "low-Froude-number  asymptotic  approximation"  to  the  surface  integral  in  formula 
(10),  we  may  obtain  the  following  low-Froude-number  asymptotic  approximation  to  the 
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-2  2 

- r -IF  (5cos9  + nsin0)sec  9 2 2 

nQ(e)  ^ F e (cos  0 - v )vdo  as  F - 0 . (17) 

' c 

In  the  low-Froude-number  limit,  the  exponential  factor  in  the  above  integral  may  be 
seen  to  be  rapidly  oscillatory,  so  that  the  major  contribution  to  this  integral  as 
F -*•  0 stems  from  the  "points  of  stationary  phase",  which  are  defined  by  the  equation 
E,'c os0  + n^sin©  = 0,  where  E,'  = d£/d a and  r\'  = di"|/do.  At  a point  of  stationary 
phase  we  then  have  tan0  = - £' /r\'  = - d£/dn  = - dx/dy,  so  that  we  may  obtain 
cos^9  = 1/(1  + tan^O)  = dy^/(dx^  + dy^)  = (dy/ds)^,  and  the  factor  cos^9  - in 
the  integral  (17)  vanishes,  since  we  have  u = n'  • i = -dy/ds.  It  may  therefore 
be  seen  that,  in  the  low-Froude-number  limit,  the  major  contributions  of  the  surface 
integral  and  of  the  line  integral  in  formula  (10)  for  the  "zeroth  approximation" 

£2^(0)  exactly  cancel  out  each  other,  which  evidently  results  in  a drastic  reduction 
in  the  value  of  the  wave  resistance  R. 

In  summary,  the  main  result  presented  in  this  section  is  the  "zeroth  approxi- 
mation" £)q(9)  to  the  "Kochin  free-wave  spectrum  function"  £2(0) . This  new  approxi- 
mation £2q,  which  is  given  by  formula  (10),  may  be  regarded  as  a generalization  of 
approximations  obtained  previously  by  Michell  in  1898,  Hogner  in  1932,  and  Maruo, 
Tuck,  and  Vossers  in  the  early  1960's;  more  precisely,  these  known  approximations 
correspond  to  particular  limiting  cases  of  the  "zeroth  approximation"  £2q,  as  it  was 
shown  explicitly.  In  particular,  the  zeroth  approximation  £2q  differs  from  the  Hogner 
"fine-ship  approximation"  £2^  given  by  formula  (11)  by  a line  integral  around  the 
ship  "waterline".  While  this  line  integral  is  small  compared  with  the  Hogner  surface 
integral  in  the  "fine-ship  limit",  that  is  as  the  angle  a between  the  "waterline"  (c) 
and  the  x axis  vanishes,  and  would  indeed  be  ignored  in  a systematic  "fine- (or  thin-) 
ship  perturbation  analysis",  the  line  integral  is  actually  equal  to  the  Hogner 
surface  integral,  at  least  to  first  order,  in  the  "low-Froude-number  limit",  that  is 
as  the  Froude  number  F vanishes,  as  it  was  just  shown.  It  may  therefore  be  seen 
that  we  have 

lim  R (F,  a = 0)  j 11m  R (F  = 0,  a)  , 

F -*■  0 U a -*•  0 U 


that  is  the  limit  a -*■  0,  F + 0 is  not  a uniform  limit.  This  interesting  property 
clearly  is  related  to,  indeed  is  at  the  origin  of,  the  "low-Froude-number  non- 
uniformity" of  the  classical  "thin-ship  perturbation  theory". 
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4.  The  first  approximation  R^ 


After  the  "zeroth  approximation"  &q(9),  Rq  — which  corresponds  to  merely 
taking  <p  as  zero  in  formula  (7)  for  the  Kochin  free-wave  spectrum  function  (2(0)  — 
an  evident  next  level  of  approximation  to  the  Kochin  spectrum  function  £2(9)  is 
obtained  by  using  the  initial  potential  given  by  formula  (2)  as  approximation 
to  the  disturbance  velocity  potential  (j>  in  the  ship  "near  field"  in  formula  (7) 
for  £2(0).  This  yields  the  "first  approximation"  £2^(0),  which  thus  is  given  by 


£2,(0)  - S2.(0)  + f <t>XE  da  +*S  [E°(a0I  + x<J>*)  - (J>IE°]uds  + 

1 0 Jh  n Jc  s t x 

+ ( f0<lnl  dxdy  + ( Eqhf  da  f°qhf  vyds  ’ <18> 

' f ' h J c 

where  equation  (10)  was  used,  and  the  notation  = 4>T  was  used  for  convenience; 

the  "nonlinear  free-surface  flux"  q is  given  by  formula  (3a),  and  the  "hull  form 
I 

flux"  q^f  is  defined  by  equation  (lb)  with  <|>  replaced  by  <|>j.  It  will  be  noted  that 

the  various  correction  terms  for  "real-fluid  effects"  qfel,  q^,  qg,  and  q^  in  formula 

(7)  were  ignored  in  formula  (18)  for  the  first  approximation  £2^(0) , although  these 

correction  terms  could  in  principle  be  included  in  this  first  approximation. 

An  interesting  alternative  expression  for  the  first  approximation  £2^(0)  may 

be  obtained  by  applying  a common  Green  identity  to  the  functions  E and  <f)  in  the 

domain  (d^)  "inside  the  ship”,  that  is  the  domain  bounded  by  the  fictitious  hull 

surface  (h)  and  the  portion  (f^)  of  the  plane  z = 0 located  inside  the  intersection 

curve  (c)  of  the  surface  (h)  with  the  plane  z = 0.  We  begin  by  deriving  a pre- 

**► 

liminary  relation  valid  for  an  arbitrary  function  ip(x)  verifying  the  Laplace  equa- 
2 

tion  7 i(i  » 0 in  the  "interior  domain"  (d^)  defined  above  and  the  Kelvin  boundary 

condition  tp  + <p  * 0 on  (f. ).  By  applying  a common  Green  identity  to  this  function 
-►  z xx  1 

<p(x)  and  the  function  E(x;0)  defined  by  equation  (7a),  which  also  verifies  the  Laplace 
2 

equation  V E m 0 in  the  "interior  domain"  (d  ) , we  may  obtain 


OJjE  - Bp  )da  - (\pE 
h n n 


Ei(i)dxdy  , 
2 


where  the  usual  notation  E = VE*n  and  tj>  = Vip»n  was  used.  From  the  assumed  condition 

n n 

ip  + ip  “0  for  x on  (f .)  and  the  relation  E + E “0,  which  may  readily  be  verified 

Z XX  1 Z XX 

from  equation  (7a),  we  have  i^E  -Bp  *Ei|j  -ipE  m(Bp  -\pE  ) , so  that  we  may  obtain 
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(\pE  - Eip  )da 

Jh  n n 


(&P„  - \pE  ) dxdy  . 

f x x x 


By  virtue  of  a well-known  Green  identity,  we  may  then  obtain 


(<pE  ~ )da  (E^»  - ipE  )dy 

n n T x x 

J c 


-f 

J c 


(E^x  - 4)Ex)yds 


in  which  we  used  the  relation  dy  = - yds,  where  u = n'*i  as  it  was  defined  previously 
in  connection  with  the  integral  equation  (1).  Now,  we  have  \ p * 'p  s*i  + \p  t«i  + v n*i 
5 ipgCT  + + \p^\>  since  a = s*i  , T = t*i  , and  V = n*i  by  definition  [the  symbols 

l,s,t,n,cr,T,  and  V have  been  defined  previously  in  connection  with  the 
integral  equation  (1) ] , so  that  we  may  finally  obtain 


iJ;Enda  +^p  [E°(o^g  + T^t)  - i[»Ex]yds 

Vi  ' r> 


tty 


■"■da  - tp  E^iJ^vyds  , 
n n 


(19) 


where  the  notation  E^  = E(x,y,z=0;9)  was  used  in  accordance  with  equation  (7b),  and 
the  superscript  i in  the  symbol  iJj  is  meant  to  clearly  indicate  that  the  normal 
derivative  4^  is  (evidently)  to  be  evaluated  on  the  "interior  side"  of  the  surface 
(h) , which  is  important  in  the  case  when  the  normal  derivative  \p  is  discontinuous 
across  the  surface  (h) . 

By  using  equation  (19) , with  the  function  ip  taken  as  the  initial  potential  <t>T 

2 1 
[which  clearly  verifies  the  Laplace  equation  V 4>T  = 0 in  (d  ) and  the  Kelvin  condi- 

II  ^ 1 

tion  4>  + 4>  = 0 on  (f  ) ] , into  formula  (18),  we  may  obtain  the  following  alter- 

Z XX  1 

native  expression  for  the  first  approximation  fi^(0) 

^(9)  » Qq(9)  + f E(i4>^  + -4)  E°(i(f»^  + q^f)vyds  + f E°q^1  dxdy  , (20) 

J U ) Jf 


where  the  integrals  over  the  fictitious  hull  surface  (h)  and  the  line  integrals 
along  the  "waterline"  (c)  have  been  grouped  together.  From  equation  (lb)  — with 
q^£  and  <p  replaced  by  q*f  and  (fi1  , respectively  — we  have 


Tif 


^n  + qhf  * ^I^h  * ° + ^ + ^ + * n * 


where  the  symbols  h^  and  hg  refer  to  the  interior  and  exterior  sides  of  the  surface 
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(h) , respectively.  However,  we  have  the  relation 

(V0x)h  * n - WI)h  • n = i • n (21) 

i e 

as  a result  of  the  fact  that  the  initial  potential  involves  a distribution  of 
sources  of  density  V = n • i on  the  surface  (h) , as  it  may  be  seen  from  formula 
(2) . We  may  then  obtain 

+ "it  • <T  + vVh  • 5 5 <«  + • 


and  formula  (20)  may  finally  be  expressed  in  the  form 


fl1(0) 


no(0)  + 


E(v  + 4>n>Hda 


4 


E°(v  + (J)1)  vyds  + 

tl  n 


E°qnldxdy 


(22) 


It  is  interesting  to  compare  formulas  (2)  and  (3)  for  the  initial  and  second 

approximations  and  <p^  to  the  disturbance  velocity  potential  (p  and  formulas  (10) 

and  (22)  for  the  zeroth  and  first  approximations  (2q  and  (2^  to  the  Kochin  free-wave 

spectrum  function  (2;  it  may  be  seen  that  formulas  (2)  and  (10)  and  formulas  (3) 

and  (22)  directly  correspond  to  each  other  with  the  substitution  G * — * E.  Numerical 

evaluation  of  the  first  approximation  (2^(9)  defined  by  formula  (22)  may  be  divided 

into  five  basic  steps,  as  follows:  (i)  evaluate  the  zeroth  approximation  (2^(9)  defined 

by  formula  (10)  for  some  "fictitious"  hull  surface  (h) , which  may  — but  need  not  — be 

taken  as  the  wetted  hull  of  the  ship  in  position  of  rest,  (ii)  evaluate  the  initial 

potential  <p^  defined  by  formula  (2),  (iii)  determine  the  sinkage  and  trim  experienced 

by  the  sh:'p  and  the  position  of  the  "real"  ship  hull  surface  (H)  corresponding  to  the 

approximation  <pT , (iv)  evaluate  the  fluid  flux  (v  + 4*1),,  = [i  + (V<t>T).J  • N across 
l n H l H 

the  ship  hull  (H)  and  the  nonlinear  free-surface  correction  flux  given  by  formula 

(3a),  and  finally  (v)  evaluate  the  three  integrals  shown  in  formula  (22).  The  first 

approximation  to  the  wave  resistance  R may  then  be  determined  by  evaluating  the 

Havelock  wave  resistance  integral  (9)  with  (2(9)  replaced  by  the  first  approximation 

£2j(9).  The  computational  task  involved  in  the  practical  implementation  of  the  above- 

described  successive  steps  admittedly  is  rather  considerable,  but  it  ought  however 

to  be  well  within  present-day  calculation  capabilities. 
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5.  The  second  approximation  R„ 


As  it  was  noted  previously,  the  zeroth  and  first  approximations  £)g(9)  and  £2^(9) 
may  be  regarded  as  the  first  two  approximations  in  a sequence  of  iterative  approxi- 
mations £2^(9),  k ^ 0,  associated  with  the  iterative  approximations  $q  = O,^,^,  ••• 
to  the  solution  (J)  of  the  integral  equation  (1)  for  the  disturbance  velocity  potential 
in  the  ship  "near  field".  The  second  approximation  (9)  corresponding  to  the  second 

approximation  <j>7  given  by  formula  (3)  will  now  be  derived.  By  replacing  <J>  by  — 

1 (2)  1 

which  will  be  denoted  by  for  convenience  — in  formula  (7)  for  the  Kochin  free- 

wave  spectrum  function  £}(0),  we  readily  obtain  the  following  expression  for  the  second 


approximation  f^O): 


n2(0)  = £2q(0)  + 4>  E 

h 


»da  +t 


[E° (a<J) (2)  + x<p  '2) ) - <J>(2)E°]uds  + 

St  X 


E°q(2)dxdy  + Eq<2)da 
L 4nl  J 4hf 


E°q£2)vyds  , 


where  formula  (10)  was  used,  the  ad-hoc  correction  terms  q^,  9S»  anc*  9W  Eor 

(2)  (2) 

"real-fluid  effects"  were  ignored  for  simplicity,  and  the  terms  q 1 and  q,  f are 

(2) 

defined  by  equations  (la)  and  (lb)  with  <j>  replaced  by  <)> 

(2) 

It  may  easily  be  seen  that  the  potential  4>  = <t>_  defined  by  formula  (3) 

2 (2)  L 

verifies  the  Laplace  equation  V (J)  = 0 in  the  domain  (d  ) inside  the  fictitious 

(2)  (2)  1 

hull  surface  (h)  and  the  Kelvin  condition  <j)^  + <J>^  = 0 on  the  portion  (f ^ ) of  the 

(2) 

plane  z = 0 located  inside  (h) , so  that  equation  (19)  — with  replaced  by  $ 
may  be  used.  Expression  (23)  for  the  second  approximation  f^O)  then  becomes 


fi2(9)  = Q0(9)  + 


P.i.(2)  < 


*)da  E°(i^2)  + q<2))vyds  + f E°q^2)dxdy  , (24) 

' c f 


where  the  integrals  over  the  fictitious  hull  surface  (h)  and  the  line  integrals  along 

i (2) 

the  "waterline"  (c)  have  been  grouped  together,  and  the  symbol  <)>  is  meant  for  the 

(2)  -*■  n 

normal  derivative  4>  = V<{>7  • n evaluated  on  the  interior  side  of  the  surface  (h) . 

n L (2)  (2) 

From  equation  (lb)  — with  q^  and  <J>  replaced  by  q^f  and  = 4>2 » respec- 

tively — we  have 
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where  the  symbols  h^  and  hg  refer  to  the  interior  and  exterior  sides  of  the  surface  (h) , 
respectively,  as  it  was  defined  in  the  previous  section.  However,  we  have  the  relation 

(^2>h  * n “ (^2^h  * n * 1 * n + (v  + 

i e n 

as  a result  of  the  fact  that  the  potential  involves  a distribution  of  sources  of 
density  v + (v  + 5 n • i + (i  4-  Vt})^)^  • N on  the  fictitious  hull  surface  (h)  , as 

it  may  be  seen  from  formulas  (3)  and  (2).  We  may  then  obtain 


1$n2)  + qhf}  * (V  + + (V  + <<)n2))H 


(25) 


(2)  (2) 

We  may  express  formula  (la)  — with  and  0 replaced  by  q and  cf) y ' — in 


the  form 


(2)  (2)  ,.(2)  ..(2), 
qni  = m - (4»z  + *xx  )z  . 0 


(2) 

where  the  term  ir  thus  is  defined  as 


.(2) 


4>z(2)  +4>^)  + <N(2)I2>x  + f^(2)*vN(2)|2 


z = - <}) 


(2) 


- tN 


(2)  | 2 


(26) 


However,  we  have  the  relation 


(Of  + ♦«>)*  = 0 = * qnl  for  (x’y)  on  (f) 


(27) 


as  a result  of  the  fact  that  the  potential  involves  a "free-surface  flux"  distri- 

X * 

bution  of  strength  q^  on  (f),  as  it  may  be  seen  from  formula  (3).  We  may  then  obtain 


I 

“ill 


+ TT 


(2) 


(28) 


By  substituting  equations  (25)  and  (28)  into  expression  (24)  for  fijO),  and  by 
using  expression  (22)  for  fi^(0),  we  may  finally  obtain  the  following  alternative  ex- 
pression for  the  second  approximation  SL(9) 


n2(0) 


fljte)  + 


E(v  + <j/2^)  da 

n n 


-t 


E°(v  + <j/2^)  vyds  + 

n n 


E°Tr^^dxdy 


(29) 


which  is  identical  in  form  to 


expression  (22)  for  the  first  approximation  fi^(0). 


■1 
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6.  Discussion 


Higher  approximations  R^,  k >_  3,  to  the  wave  resistance  R can  easily  be  defined 

in  principle.  Specifically,  it  can  be  shown  that  the  approximations  (J>  = <j>^  and 

(2)  1 

$2  = 4>  to  the  disturbance  velocity  potential  <j>  defined  by  formulas  (2)  and  (3) 

actually  correspond  to  the  first  two  approximations  in  the  sequence  of  iterative 
(k)  . 

approximations  $ defined  by  the  "zeroth  approximation" 


*“»  i 0 


(30a) 


and  the  recurrence  relation 


, (k+1) 


(xQ)  = <0 


(k) 


(x.)  + 


G(vWk))  da 

n ti 


-f 

J c 


G(v-Hj) (k) ) vyds  + 
n H 


(k) 


Git 


(k) 


dxdy  , k>0,  (30b) 


where  the  nonlinear  free-surface  correction  flux  ttv  ' is  given  by 


.(k) 


-E 


(k) 


(k) 


+ ^;'  + <|v*(k)i2>x  + >(k).v|v* 


(k)  1 2 


2 


(30c) 


The  approximations  (2q,  (2^,  and  ^ to  the  Kochin  free-wave  spectrum  function  (2(9) 
similarly  correspond  to  the  first  three  approximations  in  the  sequence  of  iterative 
approximations  (2^(0)  defined  by  the  "zeroth  approximation" 


n0(e)  = 


Evda 


- 1 E°' 


2 . 
v yds 


(31a) 


associated  with  the  "zeroth  approximation"  <f>^  = 0 to  $ , and  the  recurrence  relation 


V9)  = ^k-l(9)  + 


E(v-H|/k))Hda  -4  E°  (v+4)fk)  )TIvyds  + 

' c 


n 


H 


E°TT(k)dxdy,  k>l  . (31b) 


The  Havelock  wave-resistance  formula  (9)  — with  (2(9)  and  R replaced  by  (2^(9)  and  R^, 
respectively  — then  readily  defines  the  sequence  of  iterative  approximations  R^,  k>0, 
to  the  wave  resistance  R.  To  the  above  equations,  we  should  actually  add  equations  for 
determining  the  lift  and  moment  and  the  resulting  sinkage  and  trim  experienced  by  the 
ship.  It  may  also  be  noted  that  although  the  ad-hoc  correction  terms  for  "real-fluid 
effects"  qbl,  q^,  qg,  and  q^  in  the  integral  equation  (1)  and  formula  (7)  were  ignored 
for  simplicity  in  the  iterative  approximations  <t>^  and  (2^  defined  by  equations  (30) 
and  (31),  these  correction  terms  could  be  incorporated  in  principle.  The  iterative 
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scheme  defined  by  the  recurrence  relations  (30b)  and  (31b)  can  be  pursued  in  principle 

(assuming  convergence)  until  the  "hull  flux"  (v  + <j/k^)  = [i  + (V<t>^)  ] . jj  ancj  the 

(k)  n H H 

"nonlinear  free-surface  flux"  it  are  sufficiently  small,  that  is  until  the  boundary 

conditions  at  the  ship  hull  surface  and  at  the  free  surface  — which  may  indeed  be 

expressed  in  the  form  [i  + (V<J>)  ] • N = 0 and  tt  = 0,  respectively  — are  verified 

n 

within  desired  accuracy. 

Whereas  in  principle  the  above-defined  approximations  $ , 5^,  can  readily  be 
determined  by  using  the  recurrence  relations  (30b)  and  (31b)  and  the  Havelock  wave 
resistance  formula  (9) , the  enormous  computational  task  involved  in  the  actual  numer- 
ical evaluation  of  these  higher  approximations  drastically  restrict  the  feasibility  of 
the  approach  in  practice.  As  a matter  of  fact,  evaluation  of  the  second  approximation 
<|>2»  ^2,  R2  may  already  be  too  considerable  a task  in  practice,  and  the  first  approxi- 
mation (Jij,  R^  defined  by  formulas  (2),  (10),  (22),  and  (9)  may  well  be  the  point 

up  to  which  it  is  actually  feasible  to  pursue  the  iterative  scheme  defined  by  the 
recurrence  relations  (30b)  and  (31b).  For  practical  purposes,  the  first  approximation 
R^  may  indeed  be  regarded  as  the  main  result  of  the  present  study. 

Alternative  methods  of  refining  the  first  approximation  <J>j,  R^  that  are 
simpler  to  implement  numerically  than  the  continuation  of  the  above-described  iterative 
scheme  beyond  the  first  approximation  may  however  be  envisioned.  Two  such  methods 
may  indeed  be  mentioned  here  briefly  (as  they  will  be  described  in  detail  in  Parts  4 
and  5 of  the  present  potential  theory  of  steady  motion  of  ships) . A very  simple  modi- 
fication of  the  first  approximation  <J>  which  seems  likely  to  improve  somewhat  this 
approx  '.mat ion  is  given  by  the  "modified  initial  potential"  <t>'  = k^,  where  k is  a 
constant  which  can  be  determined  from  a straightforward  consideration  of  potential 
flow  about  a triaxial  ellipsoid,  with  main  dimensions  equal  to  that  of  the  given  ship, 
in  the  "zero-Froude-number  approximation"  (in  which  the  free  surface  is  replaced  by  a 
rigid  wall),  as  it  will  now  be  explained.  The  "zero-Froude-number  integral  equation", 
that  is  the  integral  equation  for  the  disturbance  velocity  potential,  $ say,  in  the 
"zero-Froude-number  approximation",  takes  the  form 


<H(x0)  = tfjOcg)  + f [t|)(x)  - iJ)(Xq)  ]G®(Xg,x)da(x)  , 
J h 


(32) 


where  the  function  ^(Xg)  is  the  "zero-Froude-number  initial  potential"  defined  as 


G (xg,x)v(x)da(x)  , 


(32a) 


and  G^(x0,x)  is  the  "zero-Froude-number  Green  function"  given  by  47rG^(xn,x)=-l/ |x-xn  | 
-l/|x-x^|,  with  x^  representing  the  mirror  image  of  the  point  x^  with  respect  to  the 
plane  z=0.  Clearly,  the  integral  equation  (32)  corresponds  to  the  integral  equation 
(1) , with  the  various  ad-hoc  correction  terms  for  "real-fluid  effects"  being  ignored, 
while  the  initial  approximation  i|)j  given  by  formula  (32a)  corresponds  to  the  approxi- 
mation 4)^  given  by  formula  (2).  It  is  interesting  that  in  the  case  of  a triaxial 
ellipsoid,  say  with  semiaxes  of  length  a,  b,  and  c,  the  "modified  initial  potential" 
ifl'  defined  as  \j)'  = — where  k is  some  given  function  of  a,  b,  and  c,  or  more 

precisely  of  the  ratios  b/a  and  c/b  — actually  is  the  exact  potential  ij)  . This  re- 
markable result  suggests  that  the  "modified  initial  potential"  <p'  - kup^,  where  k is 
the  above-mentioned  function  of  the  ratios  t/a  and  c/b  (that  is  the  beam/length  and 
draft/beam  ratios) , may  provide  some  improvement  of  the  original  initial  potential 
given  by  formula  (2) . The  "modified  initial  potential"  4>j,  and  the  associated  modified 
approximations  and  R^  will  be  examined  in  some  detail  in  Part  5 of  this  study.  One 
may  also  seek  to  improve  the  initial  potential  <J>^,  and  consequently  the  approximations 
and  R^ , by  taking  advantage  of  the  fact  that  the  "fictitious"  hull  surface  (h)  in 
formula  (2)  may  be  chosen  at  will,  at  least  to  a certain  extent.  An  obvious  way  of 
exploiting  this  arbitrariness  in  the  choice  of  the  surface  (h)  is  provided  by  the 
"method  of  coordinates  straining",  which  was  indeed  used  previously  for  a quite 
analogous  purpose  by  Noblesse  and  Dagan  [10].  This  modification  of  the  analysis 
developed  in  [10]  will  also  be  presented  in  Part  5. 

Besides  seeking  to  refine  the  first  approximation  cf) ^ , 0^,  R^  we  may  also  seek  to 
simplify  this  first  approximation,  particularly  from  the  point  of  view  of  numerical 
implementation.  Two  such  simplifications  may  also  be  briefly  mentioned  here.  A 
notable  simplification  of  the  above-defined  first  approximation  may  be  achieved  in 
the  case  of  ships  operating  at  fairly  low  values  of  the  Froude  number,  for  which  the 
"zero-Froude-number  initial  potential"  given  by  formula  (32a) , or  better  yet  the 
modified  potential  = kij)j,  may  be  used  as  an  approximation  to  the  potential  <J>  in  the 
ship  "near  field"  in  formula  (7)  for  the  Kochin  free-wave  spectrum  function  fi(9), 
thereby  defining  a "low-Froude-number  slender-ship  approximation".  This  approximation 
will  be  examined  in  detail  in  Part  4 of  this  study,  which  will  be  concerned  with  "low- 
Froude-number  approximations".  A different  sort  of  simplification  of  the  first 
approximation  (p^.,  can  be  achieved  by  again  exploiting  the  fact  that  the 

"fictitious"  hull  surface  (h)  may,  to  a certain  extent,  be  chosen  at  will.  Specifically, 
important  simplifications  can  be  achieved  by  assuming  the  framelines  of  the  fictitious 
hull  (h)  to  be  trapezoidal  in  shape,  as  it  will  be  shown  explicitly  in  Part  5. 
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A comparison  between  the  classical  Michell  "thin-ship  approximation"  R^  and  the 

new  "zeroth  approximation"  Rq  obtained  in  the  previous  section  is  shown  in  figure  2 

for  the  case  of  the  wedge-like  ship-bow  form  defined  by  the  equations  y = ± xtana  for 

0 <_  x <_  t and  0 >_  z >_  - ao,  and  y = ± £tana  for  t <_  x <_  °°  and  0 ^ z > - °°.  This  ship- 

bow  form  thus  is  of  infinite  draft,  and  consists  of  an  "entrance  wedge"  of  angle  2a 

2 

and  (dimensionless)  length  £(ELg/U  where  L is  dimensional)  continued  by  two  parallel 
vertical  walls  extending  to  infinity  downstream.  Specifically,  figure  2 shows  the 

function  Rn(£;a)  for  0 <_  t £ 20  and  for  a = 0° , 10°,  20°,  30°,  40°,  50°,  and  60°, 

* u *_2  ,26 
where  Rq  is  defined  as  Rq  = R^/tan  a,  with  Rq  dimensionless  (we  have  R = R'g  / pU 

where  R^  is  dimensional).  As  it  was  discussed  previously,  the  "zeroth  approximation" 

Rq  reduces  to  the  Michell  approximation  R^  in  the  "thin-ship  limit",  that  is  as 

a -*■  0.  This  may  indeed  be  verified  from  figure  2 where  the  curve  corresponding  to 

a = 0°  is  actually  identical  to  the  Michell  approximation  R^(£)  = R^/tan  a = 

Rg(£;a  = 0°),  and  both  the  symbols  0°  and  M (for  Michell)  are  indeed  attached  to  this 

curve.  It  will  be  noted  that  the  Michell  approximation  R^£,a)  is  proportional  to 

2 * 2 ™ * 
tan  a,  that  is  we  have  R^-d.a)  » RM(£)tan  a,  and  the  curves  R (■£; a)  therefore  collapse 

into  one  single  curve  R^(-£)  in  the  Michell  approximation,  while  we  do  have  a family 

of  different  curves  in  the  "zeroth  approximation"  Rq  . 

It  may  be  seen  from  figure  2 that  the  curves  R*(£;a)  generally  lie  below  the 
* * 0 

Michell  curve  ^(c)  = Rq C-c ;ot  =*  0 ) and  that  the  differences  between  the  "zeroth  appro- 
ximation" Rq  and  the  Michell  approximation  R^  are  quite  significant  even  for  such 
relatively  small  entrance  angles  as  a = 10°  and  a = 20° , and  are  in  fact  rather 
spectacular  for  larger  entrance  angles,  say  for  a _>  30°.  In  particular,  the  amplitude 
of  the  oscillations  in  the  Michell  wave-resistance  curve  are  notably  reduced  in  the 
"zeroth-approximation"  even  for  such  relatively  small  angles  as  a * 10°  and  a = 20°, 

while  for  larger  angles,  say  for  a > 30°,  these  oscillations  are  drastically  reduced. 

* 

One  may  also  observe  a notable  phase  shift  of  the  curves  Rn(£;a)  with  respect  to  the 
Michell  curve  K.(t) . This  phase  shift  is  directed  towards  the  higher  values  of  the 
"entrance  length"  t,  that  is  towards  the  lower  values  of  the  Froude  number  F = l , 
and  it  may  be  seen  to  increase  with  both  a and  t. 

The  trend  of  the  above  results  for  the  "zeroth  approximation"  Rq  is  encouraging, 
and  suggests  that  the  present  new  theory  may  help  explain,  and  remedy  to,  some  of  the 
typical  discrepancies  between  theoretical  and  experimental  wave-resistance  curves. 
Confirmation  of  this  hope  must  however  await  further  numerical  calculations,  and 
comparisons  with  experimental  results.  To  this  end,  numerical  calculations  for  the 
case  of  a parabolic  strut  are  currently  being  undertaken. 
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